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Abstract 

This is particularly apparent within the insects that show high variations in genitalic form between 
closely related species. Static allometry is one of the effective approaches for quantification of such 
variation. Despite the crucial roles of the parasitic Hymenoptera in the ecosystems, little is known 
about the sexual selection in this order, compared with other insect orders. We examined the 
allometry of different morphological traits in 35 males and 35 females of Habrobracon hebetor 
(Say, 1836) (Hymenoptera: Braconidae) from a laboratory colony reared in Mashhad city, for the first 
time. The aim was the investigation of allometric relationships of different body traits as a way for 
quantification of the natural selection impacts on the different body parts. 12 genitalic and 
non-genitalic body parts of H. hebetor males and females were photographed and measured. 
Principal Component Analysis (PCA) was used to explore the variance of the traits and two 
regression analysis methods to obtain the allometric slopes. All the non-genitalic traits in male and 
female wasps showed isometry except pterostigma width in male wasps which showed positive 
allometry. In male genitalia, two traits showed strong negative allometry and one trait showed 
isometry. Our findings showed that in this species males with an average size of genitalia were more 
successful in generating viable offspring than males with relatively smaller or larger genitalia sizes 
and this is irrespective of the overall body size. Our results showed for the first time that such 
stabilizing sexual selection might operate on genital size in the braconid wasps. 
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INTRODUCTION 

Extraordinary variation of genitalia structures in animals and its function in the natural selection 
and speciation is highly considered and discussed, lately. Sexual selection is discussed as the 
primary force driving such enormous variation (Leonard & Cérdoba-Aguilar, 2010; 
Cao et al. 2019). Static allometry, which is defined as the size relationship of a specific body part 
with the whole body in a population of a particular species, is one of the effective approaches for 
quantification of the natural selection impacts on the different body parts and their performance in 
an organism (Hosken, & Stockley, 2004; Ohno et al. 2003). This relationship is revealed by the slope 
attained in log-log regression of the size of a given body trait on an indicator of total body size 
(Eberhard, 1985). According to the selection hypotheses, relative size of different body parts is 
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examined by using an index for the overall body size. Several studies on genital allometry used a 
single somatic trait as the indicator of body size (e.g. pronotum length, elytra length and head 
width), although the use of a different body-size indicator may give rise to different results (Green, 
1999). A slope over 1.0 (positive allometry) indicates advantage of relatively larger traits in the 
larger individuals of the population. Slope around 1.0 (isometry) shows relatively constant size of 
the trait across the population. Allometric slopes lower than 1.0 (negative allometry) are the 
indicator of relatively smaller trait expression in larger individuals (Eberhard, 2009). 

Presumably, natural selection can adjust allometric slopes and consequently, affects the 
traits performance (Eberhard et al. 1998). Stabilizing selection is one of the main ways that 
selection can affects animals. In this type of selection, individuals with intermediate trait size will 
be favored in the population instead of exaggerated trait size which will be favored in the disruptive 
selection phenomena (Hosken & Stockley, 2004). 

According to the selection hypotheses, when trait size increases faster than overall body 
size, there is a positive allometric relationship, as seen in insect’s weapons. An isometric 
relationship exists when selection pressure is equal on the trait and overall body size, like most 
somatic traits in insects (Eberhard et al. 1998). If the selection pressure favors intermediate values 
of the trait regardless of body size, it causes negative allometry (Hosken & Stockley, 2004). 
Apparently, this relationship is common in male genital structures of insects and spiders (Eberhard 
et al. 1998; Eberhard, 2009). Conversely, Bolanos et al. (2014) found an exception to the previous 
findings. They interpreted isometric relationship of some genital traits in the damselflies as effect of 
sexual conflict situation that could change scaling pattern of genitalia in these insects. Also, some 
similar exceptions to the general pattern proposed by Eberhard et al. (1998) were previously found 
(Johnson, 1995; Cayetano et al. 2011). These findings showed that genital traits may evolve under 
different kinds of scaling patterns in different species. Moreover, different selective pressures 
which genitalia may confront during different sexual selection mechanisms need further 
investigations. 

Eberhard et al. (1998) proposed the one-size-fits-all hypothesis based on their study of 
allometric relationships of genital and non-genital morphological traits in many species of 
arthropods including insects. They indicated that all morphological traits in the male genitalia have 
allometric slopes significantly lower than 1.0 and lower than the slopes for non-genitalic body 
traits. Their finding means that sexual selection within the arthropod populations acts through 
stabilizing selection and males with the average genital sizes are more favored for gene transfer to 
the next generations. According to this hypothesis, female morphology is relatively invariable so 
that the consistency of male genital morphology is an adaptation to this situation. Subsequently, 
several studies on different orders of insects and spiders have confirmed their findings for male 
genital size (e.g., Schmitz et al. 2000; Iwahashi, 2001; Tatsuta et al. 2001; Eberhard, 2002; Tschinkel 
et al. 2003; Hosken et al. 2005; Mutanen & Kaitala, 2006; Rabieh et al. 2015; Cao & Hayashi, 2019). 
Eberhard et al. (1998) found that genitalia traits of the studied females have lower allometric 
slopes than the median allometric slope for non-genitalic traits but this difference was not 
significant. According to these results, they proposed that female genitalia may show the same 
patterns of variation as male genitalia. But they reported less variability for female genitalia than 
those of non-genitalic traits (Eberhard et al. 1998; Eberhard, 2009). Afterward, a few studies (e.g., 
Funke & Huber, 2005; Rabieh et al. 2015; Cao & Hayashi, 2019) confirmed their findings. Although, 
Nava-Bolafios et al. (2012) found that genital traits did not show negative allometry in two Odonata 
species, in both sexes. 

Despite the crucial roles of the parasitic Hymenoptera in the ecosystems and the high 
diversity of mating systems, little is known about sexual selection and sexual conflict in this order, 
compared with other insect orders. Although, the mating systems of this diverse taxonomic group 
are studied in-sight of how population structure influences sex allocation strategies and how to 
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encourage their attack and control of pest populations as one of the main groups of biological 
control agents in the applied research (West, 2009; Boulton et al. 2015). 

Allometric relationships of different body traits including genitalic parts are poorly studied 
in the Hymenoptera species and most allometry studies are focused on allometric relationships of 
non-genitalic traits, especially in the case of sexual dimorphisms investigations (e.g., Perrard et al. 
2012; Benitez et al. 2013; Polilov & Makarova, 2017; Quezada-Euan et al. 2019). To determine 
whether genitalic parts are coordinating with the stabilizing selection prediction of strong negative 
allometry of genitalic structures, in Braconid wasps, this paper examined the allometry of different 
genitalic and non-genitalic morphological traits in males and females of Habrobracon hebetor (Say, 
1836) (Hymenoptera: Braconidae). This species is a polyphagous and gregarious ectoparasitoid 
wasp which is widely used in integrated pest management programs. This species is mass-reared 
for control of the larvae of some important and well-known Lepidoptera pest species (Borzoui et al. 
2016). The present paper also examined CV (coefficient of variation), CV’ (degree of dispersion of 
points around the allometric line) and SEE (standard error of estimate) that are other sources of 
information about the trait's variation. The results are discussed in evolutionary frameworks. 


Material and Methods 
Insects 
A laboratory stock colony of Habrobracon hebetor was established from individuals collected from 
tomato fields during spring 2018 around Mashhad, Khorasan-e Razavi province of Iran. Rearing 
was done in the laboratory conditions at 25+1 °C, 65+5 % RH, and a photoperiod of 16:8 (L:D) on 
fifth-stage larvae of Ephestia kuehniella for at least 10 generations for obtain to adequate 
population reared on the host at these conditions (Borzoui et al. 2016). 35 adult males and 35 adult 
females of H. hebetor were selected, randomly, from the laboratory stock colony and killed in 
Ethanol 70% solution. 

E. kuehniella, were provided according to Borzoui et al. (2016) at the above-mentioned 
laboratory conditions. Examined materials were deposited in the Insect Collection of Birjand, Plant 
Protection Department, Faculty of Agriculture, University of Birjand, Iran. 


Preparation and Measurements 

Body parts of H. hebetor males and females were dissected. In each individual, the antenna, 
forewing, foreleg and hind leg were removed and cleaned to eliminate redundant tissues. Genitalia 
of both sexes were removed from the softened surrounding tissues of the abdomen and preserved 
in Ethanol 70% solution for only some minutes prior mounting. All mentioned parts were mounted 
on Canada balsam fixative. After preparation, all body parts were photographed through a 
microscope using a C-5050ZO0OM digital camera (Olympus). 15 (three genitalic and 12 non- 
genitalic) traits in males and females were examined. For prevention of errors due to body 
asymmetry, only the right body parts were examined in all traits. The 12 measured non-genitalic 
traits in both sexes were the lengths of the following body parts: forewing, pterostigma, fore femur, 
fore tibia, spine of fore tibia, mid femur, mid tibia, hind femur, hind tibia, pedicel and the width of 
forewing and pterostigma (Fig. 1). The lengths of valva and ovipositor and the width of valva in the 
females and the lengths of the gonoforcep and penisvalva and the width of penisvalva in the males, 
were the measured genitalic traits (Fig. 1). TPSDIG, version 2.16 (Rohlf, 2015) was used for all the 
measurements. To minimize the measurement error, each trait was measured by the same person 
three times. Estimated measurement error for all studied traits showed relatively repeatable 
measurements (Table 1). Therefore, the average of the three measurements was used for further 
analyses. 
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FIGURE 1. Genitalic and non-genitalic traits in male and female Habrobracon hebetor. A: Forewing 
length; B: Forewing width; C: Pterostigma length; D: Pterostigma width; E: Pedicel length; F: Fore 
femur length; G: Fore tibia length; H, Fore tibia spine length; I: Hind tibia length; J: Hind femur 
length; K: Mid femur length; L: Mid tibia length; M: Valva length; N: Valva width; O: Ovipositor 
length; P: Gonoforcep length; Q: Penisvalva length; R: Penisvalva width. 


Data analysis 

We used Principal Component Analysis (PCA), conducted separately for each sex, to obtain the best 
component explaining the variance of the traits. We selected forewing as the indicator of body size 
for allometry studies in male and female wasps. Loading all the measured traits on PC1, this trait 
showed the strongest correlation among other traits. Then, Pearson’s product-moment correlation 
coefficient was calculated for the correlation between all other morphological traits and the body 
size indicator. 

To estimate the allometric slopes, RMA and SMA regression analysis method was performed 
for log10-transformed values of the traits on log10-transformed forewing length in each sex 
according to Rabieh et al. (2015) (Table 2) 

Also, CV (coefficient of variation), CV’ (degree of dispersion of points around the allometric 
line) and SEE (standard error of estimate) that are other sources of information about the trait's 
variation were calculated (Rabieh et al. 2015). The estimated components were compared between 
genitalic vs. non-genitalic traits by Mann-Whitney U-tests within each species. All statistical data 
analyses were performed using SPSS version 19.0 (IBM Corp, 2010). 
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TABLE 1. Measurement error and coefficients of variation (CV) evaluation for all male and female 
traits in Habrobracon hebetor 


Trait Measurement CV 
Error (%) (%) 
yn-genitalic 

Forewing length 0.78 2.11 

Forewing width 1.16 2.66 

Pterostigma length 1.25 1.93 

Pterostigma width 0.99 3.23 

Fore femur length 0.27 2.19 

Fore tibia length 0.82 1.5 

Fore tibia spine length 2.13 2.31 

Mid femur length 1.17 1.64 

Mid tibia length 1.4 2.47 

Hind femur length 1.55 2.43 

Hind tibia length 0.94 1.72 

Pedicel 2.21 2.13 

Genitalic 

Valva length 1.82 2.25 

Valva width 2.01 1.3 

Ovipositor length 1.96 2.32 

Penisvalva length 2,52 1.68 

Penisvalva width 2.71 1.21 

Gonoforcep length 2.3 1.12 


TABLE 2. Allometry of evaluated morphometric traits in Habrobracon hebetor males. 


Trait r(p) ie p-Value ie a CV’ SEE 
Non-genitalic 
Forewing width 0.92(0.0008) 0.88 <0.001 0.85 0.74-1.02 1.04 0.0137 
Pterostigma length 0.88(0.0001) 0.83 <0.001 0.77 0.67-1.00 0.92 0.0159 
Pterostigma width 0.89(0.005) 1.21 <0.001 1.03 1.01-2.13 1.35 0.0217 
Fore femur length 0.85(0.0001) 0.67 <0.001 0.89 0.52-1.01 1.15 0.0145 
Fore tibia length 0.86(0.004) 0.77 <0.001 0.81 0.61-1.04 0.77 0.0167 
Fore tibia spine length 0.73(0.001) 0.78 <0.001 0.73 0.52-1.03 1.26 0.0181 
Mid femur length 0.92(0.007) 0.74 <0.001 0.65 0.63-1.22 0.64 0.0117 
Mid tibia length 0.87(0.0001) 0.75 <0.001 0.82 0.6-1.16 1.22 0.0151 
Hind femur length 0.90(0.013) 0.84 <0.001 0.79 0.7-1.31 1.06 0.0145 
Hind tibia length 0.83(0.012) 0.75 <0.001 0.69 0.57-1.12 0.96 0.0178 
Pedicel 0.71(0.005) 0.89 <0.001 0.81 0.57-1.21 1.50 0.0212 
Genitalic 
Penisvalva length 0.55(0.003) 0.2 0.006 0.18 0.06-0.34 1.04 0.0141 
Penisvalva width 0.42(0.006) 0.61 0.012 0.66 0.14-1.05 1.1 0.0157 
Gonoforcep length 0.53(0.05) 0.09 0.46 0.11 -0.16-0.34 0.94 0.0227 


-The table shows Pearson correlation coefficient (r) and its p-value, ranged major regressions 
allometric coefficient (RMA slope), the p-value obtained from permutations to test if the allometric 
coefficients differed from zero, standardized major axis allometric coefficients (SMA slope), and their 
confidence intervals (SMA slope CI), degree of dispersion of points around the allometric line (CV’) 
and standard error of estimate (SEE). 
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FIGURE 2. Principal component analysis (PCA) plots of trait loadings on PC1 vs. PC2 showing the 
covariation structure in the trait matrix of the males and females. A: Habrobracon hebetor male. B: 
H. hebetor female. FW: forewing width; FL: forewing length; F1: Fore femur length; F2: Mid femur 
length; F3: hind femur length; T1: Fore tibia length; T2: mid tibia length; T3: hind tibia length; ST1: 
Fore tibia spine length; PL; length of the pedicel; PTL: pterostigma length; PTW: pterostigma width. 
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TABLE 3. Allometry of evaluated morphometric traits in Habrobracon hebetor females. 


Trait r (p) ee p-Value ce SMA slope CI CV’ SEE 

Non-genitalic 
Forewing width 0.94(0.0002) 0.84 <0.001 0.89 0.72-1.00 0.91 0.0142 
Pterostigma length 0.89(0.0004) 0.81 <0.001 0.79 0.71-1.01 0.88 0.0161 
Pterostigma width 0.84(0.001) 0.92 <0.001 0.81 0.55-1.13 1.75 0.0188 
Fore femur length 0.79(0.0001) 0.61 <0.001 0.69 0.63-1.05 1.34 0.0154 
Fore tibia length 0.81(0.002) 0.87. <0.001 0.75 0.58-1.03 0.87 0.0173 
Fore tibia spine length 0.82(0.001) 0.6 <0.001 0.68 0.57-1.01 1.32 0.0189 
Mid femur length 0.84(0.003) 0.53. <0.001 0.52 0.43-1.00 0.88 0.0126 
Mid tibia length 0.88(0.0001) 0.68 <0.001 0.83 0.62-1.12 1.17 0.0132 
Hind femur length 0.88(0.001) 0.82 <0.001 0.77 0.7-1.25 1.15 0.0152 
Hind tibia length 0.9(0.002) 0.59 <0.001 0.64 0.58-1.02 0.74 0.0138 
Pedicel 0.81(0.004) 0.79 <0.001 0.85 0.67-1.19 1.2 0.0205 

Genitalic 
Valva length 0.48(0.05) 0.47 0.27 0.39 0.1-0.54 1.97 0.0311 
Valva width 0.69(0.01) 0.59 0.31 0.65 0.49-1.00 0.94 0.0251 
Ovipositor length 0.52(0.002) 0.62 0.01 0.8 0.53-1.01 1.98 0.0323 


- The table shows Pearson correlation coefficient (r) and its p-value, ranged major regressions allometric 
coefficient (RMA slope), the p-value obtained from permutations to test if the allometric coefficients differed 
from zero, standardized major axis allometric coefficients (SMA slope), and their confidence intervals (SMA 
slope CI), degree of dispersion of points around the allometric line (CV’) and standard error of estimate (SEE). 


DISCUSSION 
This study is a comprehensive static allometry investigation on external genitalia and non-genitalic 
traits within both sexes of a Hymenopteran species. Although, most studies on insect’s allometry 
usually focus on male genitalia traits and lay conspecific female genitalia away from their 
deduction. 

In our results, almost all non-genitalia traits in both sexes of Habrobracon hebetor showed 
isometry (greater than 1.0 in pterostigma width in males). On the other hand, male genitalia traits 
of H. hebetor showed negative allometry except for Penisvalva width which showed isometry. This 
information confirms that, in H. hebetor male genitalia and non-genitalia traits reveal coordination 
with the one-size-fits-all hypothesis proposed by Eberhard et al. (1998) and further explained in 
Eberhard (2009). This means that in a given population of this species male individuals with 
average size of genitalia are more successful in generating surviving offspring than males with 
relatively smaller or larger genitalia sizes and this is irrespective of the overall body size of them. 
Our results in terms of the negative allometry of genitalia traits vs. isometry of non-genitalic traits, 
is in coordinate with the results of other works evaluated static allometry in male genitalia of 
different insect orders (e.g., Schmitz et al. 2000; Iwahashi, 2001; Tatsuta et al. 2001; Eberhard, 
2002; Tschinkel et al. 2003; Hosken et al. 2005; Mutanen & Kaitala, 2006; Rabieh et al. 2015; Cao & 
Hayashi, 2019). Although, some studies found different results which do not support one-size-fits- 
all hypothesis, recently (e.g., Nava-Bolanos et al. 2012; Nava-Bolanos et al. 2014). These findings 
offer consideration on different functions of genitalia parts during copulation and sexual conflicts 
through different sexual selection mechanisms. For example, male genitalia parts may evolve to 
access stored sperm of females which leads to positive allometry of aedeagus in males of 
Protoneura cara (Odonata: Protoneuridae) (Nava-Bolanos et al. 2012). 

According to our results, there is no significant difference between genitalic and non- 
genitalic traits in female wasps. This finding is another confirmation for the claim that the shallow 
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static allometry slopes of male genitalic traits reflect functional roles associated with mating and 
sexual competition. 

In a given species population, larger CV of a specified trait may be due to a higher allometric 
value or a greater dispersion of points around the allometric line (CV’). The allometric slope is 
related to genetic structure of the organism which evolves under different selection pressures but 
the degree of dispersion may be related not only to genetic structure but also to other 
environmental and internal causes (Eberhard, 1998). In our results, CV’ and SEE of genitalic traits 
have values less than those for non-genitalic traits, in the male wasps. Although, their difference 
was not significant, this result means that in H. hebetor the size of male genitalia is more stable than 
that of non-genitalic parts against changes in the body size, and the difference in the degree of 
phenotypic variation between genitalic and non-genitalic morphometric feature is related to the 
difference in allometric slopes. These findings are consistent with those of previous studies which 
investigated the allometric slope and the dispersion of points around the allometric line, separately, 
for male genitalia traits and other body parts (e.g., Eberhard et al. 1998; Ohno et al. 2003; Rabieh et 
al. 2015; Cao & Hayashi, 2019). It is assumed that sexual selection in arthropods may favor a 
particular scaling relationship for male genitalia under a developmental-genetic program which 
appears as allometric slope (Eberhard et al. 1998). Our results showed that such stabilizing sexual 
selection may also operate on genital size in braconid wasps. 

Our results, presented another confirmation of the previous studies which showed that 
genitalia traits of females have lower allometric slopes than the median allometric slope for non- 
genitalic traits but this difference was not significant (e.g., Eberhard et al. 1998; Funke & Huber, 
2005; Rabieh et al. 2015; Cao & Hayashi, 2019). We found that genitalia traits in H. hebetor females 
have lower allometric slope, CV, CV’ and SEE than those for non-genitalic traits which for two last 
components, the differences were significant. Our findings showed similar allometric patterns 
among the male and female traits measured. Consequently, as male genitalia parts showed negative 
allometry and female genitalia, on the other hand displayed isometric relationship, we predict 
relatively invariability in female morphology of H. hebetor so that male wasp could adapt to this by 
being consistent in genital morphology. While our evidences do suggest the possibility of stabilizing 
selection in H. hebetor, additional experiments would be required to reinforce these claims. 

If a non-genitalic trait of a male individual attends as a secondary sexual trait during sexual 
selection mechanism, it may produce higher mating success when has a larger relative size, 
conversely, in small individuals, viability costs may limit the size of this trait, resulting in positive 
allometry (Bonduriansky, 2007). This will be much distinctive when a large number of competitive 
males exist or the trait produces resolutely more advantage in competitive interactions (Fromhage 
& Kokko, 2014). In our results, the pterostigma width in H. hebetor males showed positive 
allometry in relation to body size. Conversely, this trait in H. hebetor females showed an isometric 
relationship with body size. This suggests selection for larger pterostigma in large H. hebetor males 
and/or relatively smaller ones in small H. hebetor males. Therefore, positive allometry of the 
pterostigma width in H. hebetor may be a result of secondary sexual function. This possibility 
deserves further investigations. 


RESULTS 

The first and second components of the PCA explained 94.8% and 3.4% of the total variance in the 
males and 91.1% and 3.7% in the females, respectively. Therefore, PC1 chose as the best 
component explaining the variance of the traits for both sexes. Plots of trait loadings on PC1 vs. PC2 
showing the covariation structure in the trait matrix of the males and females are presented (Fig. 
2). Loading all the measured traits on PC1, forewing length showed the strongest correlation among 
other traits (Males: r = 0.989, p-value < 0.001; Females: r = 0.977, p-value < 0.001). All the 
measured traits were significantly related to the overall body size for male and female wasps 
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(Tables 2 and 3). In most traits, a highly significant relationship (p-value<0.01) was observed in 
both sexes (Tables 2 and 3). All the non-genitalic traits in male and female wasps showed isometric 
relationship with the body size indicator (Tables 2 and 3) except pterostigma width in male wasps 
which showed positive allometric relationship with body size (Table 2). In females, valva length 
showed strongly negative allometric relationship with body size and valva width and ovipositor 
length showed isometry (Tables 3). In male wasps, length of penisvalva and gonoforcep have strong 
negative allometry in relation to body size indicator and penisvalva width showed values which did 
not differ significantly from one (isometry) (Table 2). 

The mean slope for the measured genitalic traits of males was significantly lower than one’s 
for the non-genitalic traits (Mann-Whitney U-test, P < 0.01). In female wasps, there is no significant 
difference (Mann-Whitney U-test, P = 0.097) but the mean slope of genital traits was lower than the 
mean slope for non-genitalic traits. There is no significant difference between results of allometric 
slopes in the RMA and SMA regression methods (Mann-Whitney U-test, P < 0.01). 

Genitalic traits of male wasps showed significantly lower CV than non-genitalic traits 
(Mann-Whitney U-test, P < 0.01). Despite the lower mean value for genitalic traits, there was no 
significant difference between CV for genitalic and non-genitalic traits, in female wasps (Tables 3). 
Conversely, CV’ showed no significant difference between genitalic and non-genitalic traits in male 
wasps (Mann-Whitney U-test, P = 0.361) and it was significantly lower for non-genitalic than 
genitalic traits in female wasps. Similar situation was seen for SEE which only female wasps have 
significantly lower non-genitalic traits than genitalic traits (Tables 2 and 3). 
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